Magic Doping Fractions in High- Temperature Superconductors 
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We report hole-doping dependence of the in-plane resistivity pab in a cuprate superconductor 
La2-2,Sr2,Cu04, carefully examined using a series of high-quality single crystals. Our detailed mea- 
surements find a tendency towards charge ordering at particular rational hole doping fractions of 
1/16, 3/32, 1/8, and 3/16. This observation appears to suggest a specific form of charge order and 
is most consistent with the recent theoretical prediction of the checkerboard-type ordering of the 
Cooper pairs at rational doping fractions x = (2m -I- l)/2", with integers m and n. 
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All high-Tc cuprates contain three robust phases — 
the insulating antiferromagnetic (AF) phase, the super- 
conducting (SC) phase, and the metaUic phase — de- 
pending on the density of charge carriers introduced by 
doping. However, in some cuprate materials, there are 
also other electronic phases which compete with super- 
conductivity 1, ,2J. Determining the nature of these 
competing phases is a key focus of the current research 
in high-temperature superconductivity. One particularly 
important type of competing phase is a charge-ordered 
phase in underdoped cuprates, where the carrier density 
is smaller than the optimum level for superconductivity. 
In the underdoped regime, the mean kinetic energy of 
the carriers is low because of the small carrier density, 
and the Coulomb interaction plays an important role. 
The Coulomb interaction generally prefers some form of 
charge order, whose detailed form could be affected by 
the local antiferromagnetic exchange energy as well. One 
possibility is that the charges form one-dimensional (ID) 
stripes H, IE IE @ - Experimentally, magnetic and lat- 
tice neutron scattering on La2_a;Sra;Cu04 (LSCO) and its 
family compounds has been interpreted in terms of the 
ID stripe picture 01; although the two-dimensional 
(2D) nature of the spin system has recently been em- 
phasized in Refs. [a, 1^. Considering the presence of 
strong pairing interactions in this material, Chen et al. 
[i(t[ llll IT^ have proposed a 2D checkerboard-type or- 
dering of the hole pairs. It offers a natural explana- 
tion of the scanning tunneling microscopy (STM) results 
on BiaSraCaCuaOg+a (BSCCO) and Caa-^Na^CuOaCla 
(NCCOC) compounds, which show rotationally symmet- 
ric 4a X 4a charge ordering patterns jlji [ij, [l^ [l^ (l3| . 
The checkerboard state of the Cooper pairs has also been 
discussed in other frameworks in the recent literature 
0, [21 120|, 1^ . Furthermore, the possibility of a Wigner 
crystal of single holes has also been proposed as a com- 
peting charge ordered state at low doping 22. li^l li^ l . In 
view of the contrasting experimental results and theoret- 
ical proposals, more systematic studies of the nature of 
the charge order in the cuprates are clearly called for. 

The charge ordering tendency is expected to be par- 
ticularly pronounced near certain "magic" doping lev- 



els, where the charge modulation is commensurate with 
the underlying lattice [H HHlll H^. Motivated by 
the recent discussions on the stripe versus the checker- 
board order, we carry out a systematic study of the dop- 
ing dependence of the resistivity, in order to uncover 
the possible commensurability effects. Thanks to the 
greatly improved crystal-growth technique for LSCO us- 
ing floating-zone furnaces, single crystals of LSCO of un- 
precedented quality have recently become available 
for a very wide doping range. The cleanliness of the new- 
generation crystals has allowed, for example, to produce 
100%-untwinned single crystals |2^, which in turn led to 
finding of novel physics in this system ^27, 28, 29J. In this 
work, we systematically measure the temperature depen- 
dence of the in-plane resistivity pab in a series of high- 
quality LSCO crystals for x = 0.009 - 0.216. The raw 
data of Pab{T) are shown in Fig. 1 for all the supercon- 
ducting samples; note that the hole doping is changed in 
very small increments (typically 1%) here, which is nec- 
essary for analyzing how exactly the mobility of the holes 
changes with their density. 

Figure 2(a) shows the x dependence of the inverse mo- 
bility , which is equal to nepat, for representative 
temperatures (the hole density n is given by x/V, where 
V is the unit volume per Cu). Also, since the absolute 
values of are subject to possible geometrical- factor 
errors (which can be up to 5% in the case of our measure- 
ments), in Fig. 2(b) we show /9afc(T)/pa6(300K), which 
factors out such geometrical-factor errors. One can easily 
see that at high temperature the x dependencies of these 
variables are rather smooth and featureless, but with low- 
ering temperature they start to "oscillate" ; a peak at a; ~ 
0.13 is particularly evident. In addition, there are three 
more peaks and/or shoulders, if weaker, at a; ~ 0.06, 
0.09, and 0.18. This observation suggests that there are 
particular carrier densities where the hole motion tends 
to be hindered, which weakly enhances the resistivity at 
low temperature. Most naturally, such a behavior is in- 
dicative of a "commensurabilitv" effect associated with 
some sort of charge ordering 0, li^, 21, 22J. Remem- 
ber, in usual charge ordered systems where the Peierls 
transition is responsible, a sharp increase in resistivity is 
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FIG. 1: Temperature dependences of pab for a series of LSCO 
single crystals. The x values shown are the actual Sr contents 
measured by the inductively-coupled-plasma atomic-emission 
spectroscopy (ICP-AES). The samples are carefully annealed 
to remove excess oxygens or oxygen vacancies, so the hole 
density in the present samples are essentially equal to x. 



observed upon charge ordering |30| : in the present case, 
where the Coulomb interaction is hkely to be responsi- 
ble, the effect appears to be milder. The observed dec- 
imal numbers (0.06, 0.09, 0.13 and 0.18) suggest that 
the commensurability effect is possibly taking place at 
rational doping levels 1/16, 3/32, 1/8, and 3/16. (We 
note that there has been some prelimina.ry evidence for 
a charge ordering tendency at x = 1/16 |23.l24j'). 

Given that the inverse mobility shows an a;-dependence 
that is indicative of charge ordering, one may wonder how 
the superconducting transition temperature Tc changes 
with X] the inset of Fig. 2(a) shows the x-dependence of 
the zero-resistivity Tc in our series of samples. Besides a 
plateau-like feature for x — 0.08 - 0.12, the Tc changes 
rather smoothly without showing clear dips that can be 
associated with the magic fractions observed in the re- 
sistivity data. Hence, the putative charge order appears 
to be not particularly destructive to superconductivity; 
this is rather surprising, but is probably related to the 
STM observation [l^ that the checkerboard order shows 
up as a precursor to superconductivity. In this regard, 
it is useful to note that most of the other experiments 
concerning the checkerboard state were done below Tc 
[isl 113. ITH. Ir^ , while the existence of the magic dop- 
ing fractions is suggested in the resistivity data above 
Tc] if the charge orde ring phe nomena in cuprates involve 
the Cooper pairs 0, 0, HO, 0| and those pairs are 
formed at T > Tc, it is possible that the charge order- 
ing, observable when the superconductivity is weakened, 
is essentially of the same nature across Tc- 
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FIG. 2; (a) x-dependence of the inverse mobility fj,^^ {— 
nepab) at representative temperatures. Inset shows the 
x-dependence of zero-resistivity Tc. (b) a;-dependence of 
pab{T) / pab{300K) at T = 200, 100, and 50 K. The hole mo- 
tion tends to be hindered at low temperature at a; ~ 0.06, 
0.09, 0.13, and 0.18, which corresponds to the "magic" dop- 
ing fractions of 1/16, 3/32, 2/16, and 3/16. 



Now let us discuss the theoretical implications of our 
data. The ID stripe model predicts a particular set of 
"magic" doping fractions. The stripe model most often 
discussed in the literature involves site-centered, horizon- 
tal or vertical charge stripes separated by the AF do- 
mains at a commensurate distance of d = pa, where p 
is an integer and a is the lattice constant. (For the case 
of p = 4, see, for example. Fig. 1 of Ref. 0). The 
holes fill alternating sites on the charge stripe. In the 
stripe literature, it is commonly assumed that the hole 
doping on the stripe stays fixed, while the inter-stripe 
separation varies to accommodate different values of the 
doping level. This simple picture predicts magic doping 
fractions of a; = l/2p, with a charge unit-cell of 2a x pa. 

On the other hand, the 2D checkerboard-type order 
generally leads to a different set of magic doping frac- 
tions. Stimulated by the checkerboard orders observed 

by STM IilElElEE3> 

a global phase diagram of 
cuprate superconductors has been theoretically proposed 
and numerically analyzed 0, 0, 0|, where the zero- 
temperature phase diagram was examined in the two- 



dimensional parameter space of chemical potential ver- 
sus the ratio of the hole kinetic energy over the Coulomb 
interaction, within the framework of the 5*0(5) theory. 
Most intriguingly, this theory predicts, besides the AF 
and SC states, checkerboard-type ordering of the Cooper 
pairs at magic rational doping fractions {2m + l)/2", 
where m and n are integers A hierarchy construc- 
tion of the checkerboard states is shown in Fig. 3. The 
energetics of the checkerboard state has been studied ex- 
tensively in Ref. [l^ , both numerically and analytically. 
In general, at the magic doping fraction x = (2m-|-l)/2", 
the charge unit-cell is 2("+^'/^a x 2("+^^/-^a, pointing 
along the original Cu-0 bond direction when n is odd, 
and along the diagonal direction when n is even. This 
theory relies on mapping the original fermionic de gree s 
of freedom into effective bosonic degrees of freedom ; 
such mapping may be justified in the underdoped and 
optimally doped regimes, but fails in heavily doped sam- 
ples. Therefore, while the bosonic theory predicts all 
magic doping fractions at x = (2m -I- 1)/2"', one can only 
expect the effective theory to be valid for x < 1/4 = 25%. 
Also, it is generally expected that the charge ordering 
tendencies are stronger at higher levels of the hierarchy, 
with smaller n. 

Both the ID stripe model and the 2D checkerboard 
model adequately explain the dominant magic doping 
fraction at a; = 1/8, but they predict different sets of 
magic doping fractions, at which the system is expected 
to develop charge ordering tendencies. In this regard, 
our extensive data set on the doping dependence, if it 
indeed reflects a charge ordering tendency, contains suf- 
ficient accuracy to distinguish between the simple ID 
stripe model and the 2D checkerboard model discussed 
above. The simple stripe model predicts commensurate 
effect at magic doping fractions 1/4, 1/6, 1/8, 1/10, 1/12, 
1/14, 1/16, which should be either equally strong or vary 
monotonically in strength; therefore, the absence of any 
commensurability effects at a; = 1/6 and x = 1/10 in our 
data or any other previous experiments is puzzling in the 
simple stripe model. Although the stripe structure can 
yield a complex "devil's staircase" of commensurate dop- 
ings in nickelates [IJl ) the magic fractions suggested here 
in a cuprate superconductor would be a challenge to the 
stripe picture. On the other hand, the suggested series 
(1/16, 3/32, 1/8, and 3/16) agrees surprisingly well with 
the magic doping fractions predicted from the checker- 
board model discussed above, up to the level n = 5. At 
this level, the absence of the 1/32 fraction is understand- 
able, since the hole-pair lattice at this doping fraction 
would be very dilute and therefore disordered. More no- 
table is the absence of the 5/32 fraction, for which we 
do not have an adequate explanation at this moment. In 
passing, let us briefly mention the Wigner crystal of sin- 
gle holes, which has the charge unit-cell 2"/^a x 2"/^a at 
X = 1/2". The main difference between the pair checker- 
board and the hole checkerboard is not only the size of 
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FIG. 3: Hierarchical construction of the checkerboard-type 
ordering of the hole pairs at "magic" doping fractions (2m -|- 
l)/2", where m and n are integers. Following the construc- 
tion of Ref. [Till , the original Cu02 lattice is grouped into 
non-overlapping plaquettes, which can be represented by the 
squares on a checkerboard. A checkerboard can be alternately 
colored black and white; in our case, each black square con- 
tains four sites and no holes, while each white square contains 
four sites and two holes in the form of a Cooper pair. Such a 
state has hole doping density a; = 1/4 (a), as represented at 
the highest level of the hierarchy (e). (Electrons are denoted 
by black dots, and holes are denoted by open dots; since we 
only address the issue of charge ordering here, the spin of 
the electron is not explicitly indicated). At the next level of 
the hierarchy, consider the lattice of white squares only, and 
alternately color half of them black. Such a state has hole 
doping density a; = 1/8 (b). At one further level down in 
the hierarchy, one can either consider the lattice of the white 
squares, and alternately color half of them black, thus obtain- 
ing a state with a; = 1/16 (c), or one can consider the lattice 
of the newly colored black squares and alternately color half 
of them white, thus obtaining a state with x = 3/16 (d). This 
hierarchy construction can be obviously iterated ad infinitum, 
generating a binary tree of magic doping fractions as shown 
here in (e). 



the charge unit-cell but also their orientations, which are 
45° with respect to each other; thus, the two should be 
easily discernible in experiments at a given doping. 

Neutron scattering on the LSCO-based materials have 
provided an extensive set of data on the spin order. We 
note, however, that the nature of the spin order may 
not be directly related to the nature of the charge or- 
der reflected in the transport properties, particularly in 
the superconducting doping re gim e of LSCO where the 
spin order is mostly dynamic |32l| . If the magnetic in- 
commensurability observed by the neutron scattering is 
part of some dynamic dispersing mode |3, IS i it is 
natural that the incommensurability does not represent 
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the unit-cell of the incipient charge order. Furthermore, 
at cc = 0.02, the magnetic neutron scattering found static 
and unidirectional spin stripes that are no less than 30 
unit-cell apart (magnetic incommensurability S is 0.016) 
|34| . which would cause a large resistivity anisotropy if 
the charges are conforming to the ID spin stripes. How- 
ever, transport measurements on untwinned single crys- 
tals have found only a factor of 1.5 resistivity anisotropy 
between the "longitudinal" and "transverse" directions 
|27l |: this is rather difficult to understand without invok- 
ing some 2D character (which can, for example, be com- 
ing from a nematic stripe order [s^) in the charge system. 
In passing, we note that our transport measurements of 
LSCO have found evidence for charge self-organization 
Is^ and modest one-dimensionality [i^, E^l only in 
the lightly doped regime {x < 0.05) of this compound, 
while the present study suggests the existence of the 
checkerboard order only in the superconducting doping 
regime (x > 0.06). We also note that the measurements 
of the in-plane resistivity anisotropy become impractical 
for X > 0.06, because the structural transition temper- 
ature (below which the system becomes orthorhombic) 
comes close to or below the room temperature at these 
dopings, making it difficult to prepare untwinned sam- 
ples. 

There are only a few direct experimental observations 
of charge order in the LSCO-based materials. Neutron 
scattering on the Lai.48Ndo.4Sro.i2Cu04 (LNSCO) com- 
pound at a; = 1/8 reveals elastic charge order peaks js^ 
which can be interpreted either as orthogonally intersect- 
ing stripes on alternating planes or in the same plane. 
The latter case would also be consistent with the 2D pair 
checkerboard pattern with the charge unit-cell 4a x 4a. 
While it is fair to mention that the details of the charge 
order in LNSCO are more consistent with the ID stripe 
picture jl^, such a structure may well be a result of 
the ID modulation arising from the explicit symmetry 
breaking in the low-temperature tetragonal (LTT) phase 
of this particular material, and therefore may not be rep- 
resentative of the charge order in Nd-free LSCO. 

In view of the intriguing agreement of the present 
transport data with the hole pair checkerboard model, 
it would be desirable to systematically perform direct 
measurements of the charge order in the LSCO-based 
materials by some means. If the proposed checkerboard 
states are indeed realized in LSCO, the orientation of the 
charge unit cell should be along the Cu-0 bond direction 
at X = 1/8 (as is the case in BSCCO or NCCOC), while it 
should be along the diagonal direction near x = 1/16; it 
would be definitive if this rotation of the charge unit cell 
upon changing x is confirmed by a direct means. Also, 
since some fractions suggested in the present measure- 
ments are relatively weak, it would be highly desirable 
to carry out these systematic transport measurements 
under high magnetic fields or under high pressure, where 
one would expect the competing order to be enhanced 



and the magic doping fractions to be more pronounced. 
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